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1.  Introduction 


The  cylindrical  or  annular  HF  chemical  laser  known  as  the  Alpha  Laser  has  been  previously 
described.1  In  this  report,  we  will  concentrate  on  the  gain  generator  part  of  the  system,  which 
produces  the  F-atom  radicals  necessary  for  the  production  of  the  lasant,  HF*,  namely,  HF  with 
inverted  vibrational  and  rotational  state  distributions.  The  HF*  is  produced  by  the  chemical  reaction 

F  +  H2 - >  HF  +  H,  -  AH  =  31.6  kcal/moles  (1) 

The  gain  generator  is  the  cylindrical  vessel  containing  the  combustor  or  plenum  generating  the  F- 
atoms  from  a  hot,  chemically  reacting  mixture  at  a  steady  stagnation  temperature  approaching  2000K. 
The  hot  gases  are  exhausted  from  the  plenum  through  an  annular  nozzle  array  serving  as  the 
cylindrical  walls  of  the  plenum,  in  an  essentially  azimuthally  uniform  radial  flow.  The  nozzles  in  the 
array  are  a  converging-diverging  type  to  generate  a  supersonic  jet  exhaust  into  a  good  background 
vacuum.  The  trailing  edges  of  aerodynamic  blades  at  the  nozzle  exit  planes  inject  cold  hydrogen  into 
the  hot  stream.  The  gain  medium  is  formed  by  free  jet  mixing  with  the  hydrogen.  The  annular  gain 
medium  is  approximately  50  cm  in  radius  with  a  width  or  “thickness”  of  2.5-3. 5  cm.  An  annular 
resonator  structure,  the  HEXDARR  (High  Extraction  Efficiency  Decentered  Annular  Ring 
Resonator),  extracts  laser  power  from  this  medium  over  a  length  of  approximately  400  cm  in  a  double 
pass.  It  is  apparent  from  Eq.  (1)  that  laser  power  performance  will  be  strongly  proportional  to  the  F- 
atom  flow  out  of  the  nozzle  exit  planes,  but  it  will  depend  somewhat  on  the  plenum  stagnation 
temperature,  reaction  kinetics,  and  the  degree  of  F2  dissociation  at  that  plane.  The  plenum 
temperature  depends  sensitively  on  the  combustion  process,  as  we  will  see,  and  it  determines  to  a 
great  extent  the  F-atom  production  and  also  the  static  temperature  in  the  free  jet.  Both  gas  properties 
affect  the  small  signal  gain  coefficients  of  each  spectral  line  in  the  flow  direction.  Laser  power 
performance  will  also  be  affected  by  the  presence  of  small  amounts  of  F2  created  by  recombination  in 
the  plenum  and  in  the  nozzle  arrays  when  the  plenum  temperature  is  not  high  enough.  In  this  case, 
the  free  H-atom  generated  by  Eq.  (1)  can  react  with  F2  in  the  so-called  “hot  reaction”  of  the  chain  to 
raise  significantly  the  static  temperature  of  the  free  jet,  effectively  lowering  small  signal  gains  and 
even  shortening  the  annular  gain  width  by  enhancing  quenching  collision  kinetics.  The  H  +  F2 
reaction  has  an  exoergicity  of  91  kcal/mole,  almost  3  times  that  of  Eq.  (1).  Unfortunately,  in  contrast 
to  Eq.  (1),  little  of  this  large  heat  release  leads  to  productive  HF*  gain  for  lasing. 

Thus,  a  major  goal  of  gain  generator  design  has  been  to  minimize  the  presence  of  F2  in  the  free  jet. 

This  challenge  is  difficult  to  meet  because  space  laser  requirements  virtually  dictate  a  regeneratively 
cooled  plenum  of  low  mass  content.  Yet,  the  complex  structure  is  of  sufficient  size  that  practical 
laser  operation  occurs  in  a  transient  thermal  and  mechanical  condition,  which  is  difficult  to  model. 

The  empirical  Alpha  Laser  experience  has  been  essential  for  understanding  any  future  design 
approaches. 
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A  simple  model  for  predicting  F2  degree  of  dissociation  is  sought  to  guide  conceptual  design  or  test 
planning.  In  this  work,  we  propose  to  develop  a  simple,  spreadsheet  combustor  model  that  focuses  on 
predicting  F-atom  production  across  the  nozzle  exit  planes.  Our  objective  is  to  provide  sufficient 
insight  to  help  understand  observed  laser  performance  data  on  a  device  such  as  the  Alpha  Laser  and 
also  to  help  in  developing  combustor  operating  points  for  future  laser  devices.  Our  approach  is  the 
classical  application  of  Control  Volume  analysis  to  describe  the  combustor  process.2  The  advantage 
of  this  method  is  that  the  practitioner  defines  the  boundary  of  study  to  eliminate  much  of  the 
structural  and  reactive  flow  intricacies.  Once  the  Control  Volume  is  defined,  analysis  must  be 
developed  for  transient  nozzle  throat  changes  due  to  thermal  effects  and  their  linkages  to  combustor 
plenum  pressure  changes.  The  quantitative  development  of  this  linkage  is  the  one  unique  aspect  of 
this  report.  In  such  an  approach,  one  deliberately  sacrifices  the  spatial  and  temporal  details  that  might 
influence  an  actual  differential  fluid  element,  in  quest  of  averaged  quantities  that  will  reflect  the 
global  behavior  of  the  system  with  acceptable  fidelity. 

As  a  demonstration  of  its  application,  the  model  will  be  used  to  compute  the  F2  degree  of 
dissociation,  a,  as  a  function  of  time  for  two  tests  with  essentially  the  same  flow  rates  in  steady  state, 
namely,  HL902  and  HL909B.  A  short  model  uncertainty  discussion  will  follow. 
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2.  Model  Logic 


The  logic  for  building  the  thermal  energy  flux  balance  for  the  model  is  as  follows: 

1 .  Estimate  thermal  energy  flowing  in  with  reagents  during  interval  At. 

2.  Estimate  thermochemical  energy  flux  released  from  reagents,  assuming  complete 
combustion. 

3.  Estimate  net  heat  flux  loss  from  the  reagents  or  products  during  At  interval 

4.  Estimate  heat  loss  from  combustor  gases  into  the  rings  (and  end  domes)  from  temperature 
growth  in  the  rings. 

a.  Determine  growth  in  combustor  pressure  Po  (-Pc)  during  the  interval. 

b.  Relate  Pc  growth  to  nozzle  throat  area  A*  closure.  (1-D  gasdynamics) 

c.  Relate  temperature  growth  in  ring  to  A*  decrease,  (linear  coefficient  of  thermal 
expansion) 

d.  Compute  net  heat  flux  loss  from  gases  to  vessel,  given  vessel  temperature  growth. 

e.  Estimate  heat  loss  into  combustor  core  (found  negligible). 

5.  From  NET  heat  flux  from  reagents  and  product  gases  (items  1-3),  compute  “stagnation” 
temperature  Te  of  the  combustor  plenum,  (i.e.,  no  exhaust  gasdynamics). 

6.  Compute  F-atom  density  and  F2  degree  of  dissociation,  a,  given  Te. 

7.  Slightly  corrected  for  further  recombination  through  the  nozzle  structure  to  exit  plane 
(Appendix  C) 

The  main  innovation  in  this  heat  flux  model  is  to  use  the  mechanical  changes  of  the  combustor  vessel 
as  indicated  by  gasdynamics  to  deduce  vessel  temperature  growth  and,  thereby,  the  net  heat  flux  loss 
from  the  product  gas,  a  quantity  difficult  to  measure  directly.  This  vessel  temperature  growth 
manifests  the  net  heat  flux  to  the  vessel,  the  difference  between  heat  loss  from  the  product  gases,  and 
the  regenerative  cooling  of  the  vessel  by  reagents.  Since  Items  1  and  2  are  well  known,  we  will 
concentrate  on  the  development  of  Item  3  in  the  following  sections. 

It  turns  out  that  this  model  is  applicable  provided  there  are  no  mechanical  changes  in  the  vessel  due 
solely  to  Pc  pressure  changes.  We  will  discover  later  that  this  assumption  approaches  validity  beyond 
t  a  8  s  in  the  cases  discussed  here  when  the  flows  reach  a  true  steady  state  and  combustion  is  well 
established.  (See  Figure  3  and  discussion.)  An  upgraded  version  of  this  model,  which  accommodates 
the  detailed,  fast  transient  behavior  of  Pc  and  the  role  of  the  cold  inner  core  of  the  vessel,  is  reported 
in  References  9  and  10. 
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3.  Control  Volume 


We  assume  in  the. time  increment  modeled  (typically  1  s)  that  the  flow  is  “steady.”  The  combustor 
gas  properties  such  as  molecular  weight,  W,  specific  heat,  cp,  mole  fractions  of  species,  Xi?  and  ratio 
of  specific  heats,  y,  are  definable  uniform  quantities  within  a  Control  Volume  space.  These 
parameters  for  all  species  can  be  evaluated  from  the  given  reagent  flowrates  and  gas  properties  and 
those  same  parameters  deduced  for  the  product  species.  Completed  combustion  and  complete 
consumption  of  the  oxidizer  to  products  is  assumed.  There  may  be  combustor  gas  property  changes 
in  succeeding  time  intervals  due  to  slow  changes  in  the  flow  rates  of  the  reagents  and  the  observed 
change  in  combustor  pressure. 

The  choice  of  combustor  Control  Volume  is  illustrated  in  Figure  1.  We  attempt  to  choose  the  volume 
to  avoid  the  need  for  specific  three-dimensional  details  of  the  combustor  vessel  structure,  including 
the  the  heat  flux  into  reagents  during  regenerative  cooling  of  the  vessel  and  the  complex  details  of  the 
actually  mixing  of  reagents  at  an  injector  system  before  combustion. 

The  Control  Volume  suggests  the  balance  equation  for  heat  flux 


Qdotout  Qdotinflow  +  Qdotjhe^Q^gjjjj^j  Qdotgas  |oss, 


(2) 


Figure  1.  Illustration  of  the  Control  Volume  chosen.  The  Control  Volume  is 
traced  out  by  the  dotted  boundary,  which  intercepts  the  incoming 
pipes  at  the  cryo-heat  exchangers  and  the  outflow  to  the  nozzle 
throats. 
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where 


Qdotinflow  =  2j  [  mdotiCp-iTin-i] 


(3) 


With  mdot,  the  mass  flow  rate,  cp(Tin),  the  specific  heat,  and  Tjn,  the  temperature  at  inflow  of  the  ith 
reagent.  We  develop  Qdottherraochemical  and  Qdotloss  in  successive  sections  below. 
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4.  Thermochemistry 


The  combustion  process  begins  with  the  early  ignition  by  a  burst  of  the  hypergolic  mixture  of  fluorine 
and  deuterium.  The  bulk  chemicals  flowing  into  the  plenum  are  nitrogen  trifluoride,  NF3,  and 
cryogenically  cooled  (to  a  nominal  80K)  deuterium/helium  mixture.  In  the  combustion  process,  NF3 
readily  dissociates  at  elevated  temperature  to  produce  F-atoms.  A  portion  of  these  atoms  reacts  with 
D2  to  keep  the  combustor  furnace  going.  At  the  plenum  temperatures  forecasted  for  the  given  Alpha 
Laser  flowrates,  the  following  chemical  equation  governs  the  process: 

aD2  +  bNF3  +  cHe - »  2aDF  +  (l/2)bN2  +  a  (3b-2a)F  +  (l-a)(3b-2a)F2  +  cHe,  (4) 

where  coefficients  a,  b,  and  c  are  mole  fractions  based  on  incoming  flow  rates,  and  a  is  the  degree  of 
dissociation  of  molecular  fluorine. 

The  net  heat  of  reaction,  following  Eq.  (4),  is 


AHnet  —  —  aAHp2  —  bAH^^-p^  —  cAHjje  +  2aAHj}p  +  0.5bAHpj2  +  oc(3b  —  2a)AHp 

+0.5(1  -  a)(3b  -  2a)AHp2  +  cAHHe  (5) 

Of  course,  the  heat  of  formation  for  He  is  moot,  and  those  for  D2  F2  and  N2  are  zero.  The  heats  of 
formation  are  taken  at  the  condition  of  the  gas  crossing  the  Control  Volume.  AH  for  NF3  is  taken  at 
300K,  and  those  for  the  combustion  products,  at  a  themochemical  temperature  to  be  defined  below. 

A  consistent  solution  for  the  combustor  will,  therefore,  require  iteration,  in  which  gas  properties  are 
adjusted  as  we  home  in  to  a  temperature. 

We  determine 

Qdott^ermochemical  “  nc^in  #  ^^nef 

Where  ndotin  is  the  total  molar  flowrate  of  reagents.  We  also  discovered  in  the  course  of  building  this 
model  the  high  sensitivity  of  the  computed  results  to  very  accurate  thermodynamic  data  of  the  various 
species.  We,  therefore,  resorted  to  look-up  tables  based  on  the  established  JANAF  database3,4,5  and 
to  linear  interpolation  between  the  100K  temperature  intervals  of  the  database.  (Ref.  5  is  probably  the 
most  useful  for  providing  digital  data  and  functional  extrapolations.) 
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5.  Heat  Flux  Loss  Term 


The  heat  flux  loss  from  the  gases  is  given  by 

QdotioSS  =  Q^^rings  —  Q^°^struts  Q^°^domes  Q^^stagnant  core  »  ^ 

with 


Qdotrjngs  —  [Qdotconvective  Q^°Wliative  Qdotregen] 

Qdotdomes  will  become  a  small  multiplicative  correction  factor  on  Qdotfings  in  Eq.  16,  while  Qdotstruts 
and  Qdotstagnam  core  can  be  shown  to  be  negligible  heat  conductors.  Thus,  evaluation  of  Qdotloss  cen¬ 
ters  on  a  good  estimate  of  Qdotrings,  the  net  heat  flux  gained  by  the  rings. 

In  the  time  interval  At  of  assumed  “steady”  flow,  the  total  mass  flow  rates,  individual  mass  flow  rates, 
the  specific  heats,  and  the  ratio  of  specific  heats  are  all  assumed  constant.  The  construction  of  this 
model  really  assumes  that  input  mass  flow  rates  and,  therefore,  the  output  total  mass  flow  rate  remain 
constant  in  succeeding  At  intervals.  Furthermore,  the  model  assumes  that  in  steady  state,  the  vessel 
accrues  no  molecules;  namely,  the  inflow  of  reagents  equals  the  outflow  of  chemical  products.  This 
condition  is  only  approached  after  5  to  7  s  following  the  ignition  of  the  combustor.  (Lasing  begins 
about  5.5  s  after  ignition).  The  reason  is  that  key  reagent  flows  of  NF3  and  D2  must  build  up  and  are 
observed  to  asymptotically  approach  a  true  steady  state  in  an  exponential  manner.  Beyond  this  point 
in  the  actual  laser,  increases  in  pressure,  P0,  reflect  mainly  the  result  of  thermal  effects  in  the  gas  and 
the  vessel.  Then  from  the  choked  flow  equation  of  one-dimensional  gas  dynamics7  the  combustor 
pressure,  P0,  changes  because  of  two  effects:  the  square-root  dependence  on  the  stagnation 
temperature,  T0,  and  more  importantly,  the  inverse  dependence  on  the  throat  area  A*. 


P0  =  /A  * ,  where  constant  C  =  mdot 


(8) 


where  mdot  is  the  total  mass  flow  rate  out  of  the  vessel,  y  is  the  ratio  of  specific  heats,  W  is  the 
molecular  weight  after  combustion,  and  R0  is  the  universal  gas  constant.  Properties  W,  g,  and  cp 
require  calculations  combining  the  individual  properties  taken  at  the  gas  stagnation  temperature,  Te, 
of  the  product  species  according  to  their  mole  fractions.  For  example: 

W  =  [cWHe  +  0.5bWN2  +  2aWDF  +  a(3b  -  2a)WF  +  0.5(1  -  a)(3b  -  2a)WF2]ndotin/ndotou„  (9) 
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assuming  full  dissociation  into  F.  Because  the  consistent  combustor  solution  requires  solving  for  T0, 
an  iteration  on  properties  must  once  again  be  performed  to  determine  T0.  In  differential  fractional 
terms  for  rate  of  change  in  the  parameters,  Eq.  (8)  becomes 


dPo/Po  =  -dA*/A*  +  0.5dTo/To. 


(10) 


From  Eq.  (10),  one  sees  a  direct  relationship  between  a  thermal  gasdynamic  parameter,  P0,  and  a 
structural  or  mechanical  parameter,  A*,  the  throat  area,  which  will  be  affected  by  thermal  changes. 
We  develop  a  second  equation  below  to  eliminate  the  fractional  change  in  T0  (for  the  moment)  such 
that  a  fractional  change  in  P0  defines  a  fractional  change  in  A*,  and,  therefore,  the  thermal  condition 
of  the  rings  in  the  vessel  per  unit  time,  from  which  Qdotrings  can  be  derived. 

In  the  Alpha  Laser,  A*  represents  25  primary  nozzle  throats  plus  two  “half’  nozzles  on  the  ends  of 
the  ring  stack.  From  Figure  2,  we  model  one  of  the  throats  as  e  •  c,  where  c  is  the  circumference  of 
the  ring,  r  is  the  radius,  and  e  is  the  slit  width.  Due  to  thermal  effects,  both  may  change  dimension- 
ally  with  time: 


A*  =  ne  •  c  =  ns  •  27tr,  n  =  26  throats. 


(ID 


This  analysis  is  helped  by  the  fact  that  there  is  absolutely  no  thermal  expansion  allowed  in  the  direc¬ 
tion  of  the  gain  generator  cylinder’s  axis,  the  z-direction,  due  to  clamping  referenced  to  massive  invar 
rods  of  the  Alpha  Laser  structure,  as  shown  in  Figure  2b. 


ring  i+1 


Figure  2a.  Model  of  a  nozzle  throat  with  length  c  =  2nr. 


I 


I  C  =  2nr 

Figure  2b.  Concept  of  gain  generator  with  clamped  rings. 
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Due  to  the  clamping  condition,  the  change  in  area  of  A*  depends  only  on  the  thermal  expansion  of  the 
half  section  of  the  ring  cross  section  bordering  a  particular  nozzle,  and  if  one  assumes  all  the  expan¬ 
sions  are  identical  for  each  of  the  26  nozzles,  dA*  =  n  •  2rcr  •  de  +  n  •  £  •  2rc  dr. 


Between  the  rigid  centerlines  in  Figure  2a,  any  shrinkage  in  £  is  due  to  thermal  expansion  of  the  two 
material  half-lengths  Ls.  (Note  Lj  is  not  the  geometric  half  length):  d£  =  -2dL;,  and  it  follows  that 


dA*/A*  =  -2  [(dLj/Lj)  •  (L/e  )]  +  dr/r. 


(12) 


It  can  be  shown  that  L/e,  taken  constant  in  our  calculations,  is  a  very  slowly  varying  function  of  time 
t,  on  the  order  of  2%  increase  during  a  typical  test.  (See  Appendix  B:  Calculation  of  L/e.)  The 
thermal  linkage  is  through  the  coefficient  of  thermal  expansion,  CTE,  to  find  an  average  temperature 
increase  in  combustor  vessel  rings: 


dL/Lj  =  CTE  •  ATvessel.  (13) 

Using  the  perfect  gas  law  to  describe  the  combustor  gases  yields  another  relationship  for  dP0/P0. 

With  this  additional  relationship,  T0  and  r  can  be  eliminated  with  Eqs.  (10)  and  (12)  to  yield  the 
sought  relationship,  deduction  of  the  temperature  change  in  the  material  based  on  merely  the  change 
in  the  combustor  plenum  pressure: 


AT-GGA  =  (AP0/P0)  +  [4  •  CTE  •  (L/e)]  (14) 

The  estimation  of  the  heat  transfer  out  of  the  hot  combustor  gases  to  the  cylindrical  walls  of  the  vessel 
depends  on  accurate  knowledge  of  the  mass  of  a  gain  generator  ring.  (See  Appendix  A:  Mass  Data.) 
Assuming  the  convective  heat  transfer  process  at  the  boundary  of  the  end  caps  is  the  same,  one  also 
needs  a  knowledge  of  the  mass  of  the  end  caps  (i.e.,  the  domes)  of  the  vessel.  The  ring  mass  is  self 
cooled  or  regeneratively  cooled  by  the  incoming  cold  reagents.  In  the  Alpha  Laser,  the  end  caps  are 
water  cooled,  providing  a  slight  enhancement  of  heat  transfer  out  the  end  domes. 

There  is  a  detail  in  handling  the  end  domes  heat  transfer.  As  noted  in  Appendix  A,  the  end-domes 
mass  has  been  taken  as  15%  of  the  mass  of  the  rings.  The  idea  that  the  end  domes  bound  both  an 
active  annular  convective  flow  region  and  a  stagnating  core  has  been  included.  Because  the  domes 
are  water  cooled,  an  additional  multiplicative  correction,  f,  to  Eq.  16  is 


f=  {1+0.15  •  [(TAi - Twater)/(Te - Twater)] }/l.l5,  (15) 

with  TA1  as  the  ring  temperature  at  the  given  time  interval,  At.  Generally,  Twater  equals  300K. 

The  heat  transfer  rate  Qdot  (watts)  to  all  parts  of  the  aluminum  vessel  then  is  given  globally  by 
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QdotA]  =  MA1  •  Cp  _  A1  •  ATa1  •  f/At, 


(16) 


with  the  vessel  mass  being  given  by  MA1,  the  specific  heat  of  aluminum  Cp  _  A1,  ATa1  by  Eq.  (13),  f 
by  Eq.  (14),  and  the  time  interval,  At,  of  heat  transfer  for  the  AP0  observed,  namely,  1.0  s.  The  Con¬ 
trol  Volume  approach  combined  with  ideas  relating  to  Eq.  (7)  link  vessel  material  thermal  growth  to 
the  combustor  gas  loss: 

QdotA1  =  Qdotgasloss  (17) 
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6.  Calculation  of  the  Stagnation  Temperature  Te 


Finally,  by  the  heat  flux  balance  of  Eq.  (2),  gas  stagnation  temperature,  Te,  in  the  combustor  is 

Te  =  (Qdotjnflow  Qdotthcfnochemical  —  Q*^°^gas  loss)  ^  (mdot  X  Cp_gas),  (18) 

where  mdot  is  the  measured  mass  flow  rate  of  products  from  the  combustor  and  cp_gas  is  the  specific 
heat  at  constant  pressure  of  products  computed  according  to  accurate  thermochemical  data  and  mole 
fractions  defined  by  products  (e.g.  Eq.  9).  Gas  properties  are  computed  at  Te.  Qdotinf]ow  is  the  heat 
flux  into  the  Control  Volume  due  to  heat  capacity  of  reagents.  Qdotthermochemical  is  the  heat  flux 
added  to  the  flow  as  products  after  complete  combustion.  Because  the  gas  properties  of  products  are 
computed  at  Te,  iterations  are  also  needed  here  to  arrive  at  a  consistent  solution  for  Te,  including 
combustor  degree  of  dissociation  effects. 

In  fact,  it  should  be  realized  that  this  stagnation  temperature  is  really  a  mythical  “average”  exactly 
defined  by  Eq.  (18)  within  the  context  of  one-dimensional  gasdynamics  analysis.  The  real  gas  condi¬ 
tion  is  non-uniform.  A  differential  fluid  volume  may  experience  a  temperature  up  to  a  hypothetical 
“maximum”  combustor  temperature  (sometimes  called  T;)  that  could  be  computed  without  the  “gas 
loss”  heat  flux  term,  but  adding  the  heat  flux  of  the  reagents  used  in  regenerative  cooling.  In  such  a 
simple  model,  we  view  Te  as  an  average  “one-dimensional”  stagnation  temperature  of  the  gas  flowed 
through  the  nozzles  into  the  laser  cavity. 
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7.  Estimation  of  F-Atom  Production 


For  the  given  time  interval,  At,  once  Te  is  found,  the  estimation  of  F-atoms  and  molecular  F2  in  the 
combustor  follows  conventional  equilibrium  chemistry  calculations  with  Te  as  the  defined  flow 
“equilibrium”  temperature.  Inside  the  combustor,  it  is  assumed  that  flow  effects  are  negligible  and 
there  is  sufficient  dwell  time  for  the  mixture  to  equilibrate  chemically.  It  has  to  be  confirmed  that 
there  are  no  major  bulk  fluid  dynamical  recombination  effects  to  generate  F2  through  the  nozzles. 
Indeed,  there  is  a  minor  correction  of  flows  from  the  combustor  to  the  nozzle  throat  because  fluid 
particle  times  are  sufficiently  long  to  allow  some  recombination.  However,  from  the  throat  to  the 
exit,  the  flow  becomes  on  the  average  hypersonic,  and  any  effects  can  be  shown  to  be  negligible 
(Appendix  C).  It  can  be  shown  that  loss  of  F-atoms  by  wall  and  boundary  layer  recombination  effects 
is  negligible  for  this  nozzle  design,  based  on  wall  data  discussed  by  N.  Cohen  et  al.6  Thus,  the  prin¬ 
cipal  loss  of  F-atoms  is  presumed  to  occur  by  three-body  bulk  recombination  in  the  presence  of  dilu¬ 
ents  and  other  exhaust  products. 

The  dissociation  of  F2  into  F-atoms  in  this  model  is  given  by 

kf  =  5  x  1013  a  exp[-35,000/RoTe],  (cm3-mole_1-s_1  ),  (19) 


a  =  1  for  all  collision  partners  except  He,  a  =  2,  F2  itself,  a  =  2.7,  and  F-atoms,  c  =  10.  The  specific 
reference  is  N.  Cohen  and  J.  F.  Bott.7  Instead  of  the  more  commonly  JANAF-accepted  bond  strength 
of  38  kcal-mole-1>  we  have  used  35  kcal-  mole-1  as  recommended  by  Cohen.  The  backward  rate 
coefficient  for  this  model  is  fitted  to  the  early  data  summary6 

kb  -  (2.6  x  1016)/Te105,  (cm6-mole-2  s-1)  (20) 

and  the  equilibrium  constant  K  is  simply  kf/kb.  The  maximum  number  of  F-atoms  (i.e.,  the  F-atom 
density  [F0])  available  in  the  combustor  can  be  computed  following  Eq.  (4)  with  unity  degree  of  dis¬ 
sociation  and  for  very  high  T  ,  given  the  gas  properties  of  the  gas.  For  the  1-s  interval  of  apparent 
steady  state,  this  is  equivalent  to  the  flowrate  of  F  atoms  that  will  leave  the  plenum  as  well.  However, 
the  actual  number  of  F-atoms  in  the  combustor  that  makes  up  the  desired  F-atom  flow  rate  depends  on 
T  .  The  solution  for  the  actual  F-atom  density  [F]  is  quadratic,  for  the  positive  branch 

[F]  =  0.25  •  K{-l  +  ^/l  +  8F0  /K}  (21) 

It  follows  that  in  the  given  local  fluid  particle,  density  [F2]  and  degree  of  dissociation,  a,  are 


[F2]  =  0.5  { [  F0]  -  [F] } 


(22) 
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a  =  {l  +  2  •  ([F2]/[F])}-1 


(23) 


The  admittedly  large  uncertainty  in  the  activation  energy  of  Eq.  (19)  and  the  empirical  curve  fit  in  Eq 
(20)  mean  the  uncertainties  in  the  estimations  of  F-atom  and  F2  densities  will  remain  high  except  in 
the  vicinity  of  unity  a.  Model  uncertainties  will  be  discussed  further  below.  The  efficient  HF  con¬ 
tinuous-wave  chemical  laser  usually  requires  the  suppression  of  F2  species.  This  is  achieved  by 
ensuring  sufficiently  high  Te,  in  the  combustor.  For  a  space  system,  a  requirement  for  high  Te  must 
be  traded  against  the  mass  increase  in  reagents  for  producing  the  hotter  gain  generator. 


16 


8.  Model  Results  and  Discussion 


The  model  is  demonstrated  using  recent  Alpha  Laser  input  test  data  from  Tests  HL9021  and 
HL909B8.  We  caution  at  the  onset  that  the  model  prediction  of  degree  of  dissociation  a  or  of  F-atom 
flow  rate  will  remain  unverified,  for  there  are  no  independent  direct  experimental  measurements. 
Nevertheless,  certain  observed  behaviors  of  laser  power  as  a  function  of  time  in  HL902  and  HL909B 
are  consistent  with  model  findings.  Table  1  summarizes  the  flow  condition  for  HL909B  when  the 
reagents  have  reached  a  near  steady  state  at  the  onset  of  lasing  (t  =  8  s)  In  reality,  both  NF3  and  D2 
increase  1  %  or  less  during  the  time  of  lasing.  Flow  conditions  are  similar  for  HL902. 

The  key  test  observable  is  the  measured  combustor  plenum  pressure  as  a  function  of  time  from  zero 
seconds  to  reagents  shutoff.  In  the  Alpha  Laser  test  series,  data  are  taken  at  0.25-s  intervals.  This 
work  will  use  data  only  at  each  1.0-s  interval  for  most  of  our  cases  because  the  rate  of  change  of  the 
pressure  in  this  interval  is  sufficiently  small,  typically  5%  or  less.  The  estimated  errors  on  this  meas¬ 
urement  are  3%  or  less. 

From  the  HL909B  combustor  pressure  Pc  as  a  function  of  time,  one  can  deduce  the  apparent  heat  loss 
to  the  vessel  for  each  1-s  block  of  the  study.  From  the  computed  heat  fluxes  for  a  given  t,  At  block 
(Eqs.  3,  6,  16,  17),  one  can  deduce  the  probable  Te  based  on  Eq.  (18).  This  parameter  is  found  to 
have  a  considerable  range  as  a  function  of  time,  from  1423K  to  1563K  in  this  computation  during 
lasing,  from  t  =  9stot=15s.  This  spread  is  reflected  in  the  computations  for  degree  of  dissociation 
a,  as  exhibited  in  Figure  3.  (diamonds) 


Table  1 .  HL909B  Combustor  Flow  Conditions  Used  in  the  Model  (t  =  8  s) 


Parameter 

Comments 

NF3  Flowrate,  mol/s 

112.9 

Entry  at  300K 

D2  Flowrate,  mol/s 

111.5 

Entry  at  80K 

He  Flowrate,  mol/s 

573.5 

Entry  at  1 00K,  taken  as  80K 

Reagents  Molecular  Weight,  g/mol 

13.50 

Reagents  Specific  Heat  Cp,  J/g-K 

1.928 

Reagents  Sp.  Heat  Ratio,  y 

1.470 

Combustion  Products 

Complete  reaction  assumed 

Product:  max-F  (or  F0),  mole/s 

117.2 

No  recombination  or  loss 

Product:  N2,  mole/s 

56.44 

Product:  DF,  mole/s 

221.4 

Product  Molecular  Weight,  g/mol 

11.11 

Product  Specific  Heat,  cp,  J/g-K 

2.082 

MdotF0/A,  kg/m2 

0.330 

Product  Sp.  Heat  Ratio,  y 

1.561 
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Calculation  of  HL909B  Degree  of 
Dissociation 


Figure  3.  Model  estimation  of  degree  of  dissociation  for  test  HL909B. 

The  degree  of  dissociation  or  a  curve  labeled  “Model  ( lasing),”  from  9  to  15  s,  exhibits  an  approach 
to  “steady  state”  for  times  greater  than  10  s,  with  a  approaching  0.90.  Here,  the  Pc  and  reagent  flow 
rate  data  are  taken  from  test  datal  ,8  from  HL909B.  Since  data  used  are  in  1-s  intervals  of  time,  high- 
frequency  temporal  details  are  disregarded  in  this  model. 

Comparison  of  HL902  degree  of  dissociation  results  can  be  made  with  those  of  HL909B.  The  model 
can  be  run  with  HL902  pressure  and  flow  conditions,  which  are  very  similar.  HL902  a’s  (grey 
triangles,  labeled  “HL902  Model”)  are  slightly  lower  than  those  of  HL909B,  but  within  5%.  The 
good  agreement  is  confirmation  of  the  contention  that  test  reagent  flow  rates  have  been  the  same  for 
at  least  10  tests  in  the  series,  spanning  HL902  to  HL909B.  Comparison  can  also  be  made  with  a 
contractor  model  (squares,  labeled  “TRW  HL909B”)8  based  on  a  similar  heat  flux  balance  approach, 
which  uses  data  from  embedded  thermocouples  in  the  rings.  The  TRW  model  seems  to  predict  higher 
a  at  earlier  times  averaged  over  a  much  more  steady  state,  possibly  an  indication  of  insensitivity  of 
the  thermocouples  to  more  rapid  changes  in  gas  flow  conditions.  The  thermocouples  are  highly 
sensitive  to  vessel  thermal  non-uniformities  and  the  quality  of  their  placements.  The  trend  toward 
agreement  in  degree  of  dissociation  (within  10%)  of  two  model  approaches  at  late  times  past  10  s. 
suggests  progress  toward  prediction  of  this  elusive  key  parameter. 

We  have  included  a  fourth  computation  (Model  P-3)9,10,  which  attempts  to  synthesize  the  transient 
degree  of  dissociation  in  HL909B  toward  steady  state  in  the  case  of  both  pressure  and  thermal 
transient  changes.  As  we  suggested  earlier,  both  the  current  and  the  advanced  models  converge  to  the 
same  dissociation  results  as  flows  reach  a  steady  state,  at  t  >  8-10  s. 
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The  absence  of  direct  measurements  for  F-atoms  or  a  at  the  nozzle  array  exit  plane  makes  model 
validation  difficult.  We  must  resort  to  plausibility  discussions  that  demonstrate  the  model’s  consis¬ 
tency  with  other  observed  data.  We  examined  three  types  of  data: 

1 .  Typical  thermocouple  data  for  the  rings 

2.  Spectra  behavior  as  a  function  of  time  beyond  1 1  s 

3.  P(t),  outcoupled  power  as  a  function  of  time,  seen  to  be  slowly  increasing 

In  all  cases,  the  model  seems  to  provide  insight  into  data  variations  or  demonstrate  consistency. 

In  the  Alpha  Laser,  thermocouples  were  mounted  manually  in  various  gain  generator  ring  positions. 
Although  the  absolute  values  had  considerable  scatter,  they  seemed  to  exhibit  time  constants  longer 
than  a  fraction  of  a  second,  and  a  common  trend  appeared  to  be  an  average  temperature  increase  of 
15K/s  over  most  of  the  running  time.  The  model  prediction  from  Eq.  (14)  is  shown  in  Figure  4, 
together  with  a  comparison  to  a  15K/s  slope  line.  It  can  be  seen  that  there  is  qualitative,  even  quanti¬ 
tative  agreement  over  a  wide  span  of  time. 

The  observation  of  spectra  at  times  greater  than  approximately  11  s  in  a  lasing  test  can  be  seen  from 
the  data  of  relative  laser  power  as  a  function  of  time,  as  in  Figure  5.  The  framed  region  corresponds 
to  a  lasing  time  region  when  a  has  not  reached  a  steady  state.  For  times  greater  than  around  10.5  s,  it 
can  be  seen  that  most  spectral  lines  are  steadier  as  a  function  of  time. 

The  best  demonstration  of  an  approach  to  steadiness  is  P2(6)  and  P2(7)  while  PI (7)  declines  slightly. 
Steady  spectral  signals  are  an  indication  that  laser  power  has  attained  a  steady  condition.  In  turn,  this 
implies  proportionally  that  the  F-atom  flowrate  and  therefore  a  has  become  constant  for  longer  times, 
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Figure  4.  Growth  in  ring  temperature  computed  by  the  model.  The 
line  delineates  1 5K/s. 
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Figure  5.  Outcoupled  laser  power  spectra  for  times  greater  than  11s. 

HL909B. 

as  forecast  by  our  model.  The  increases  as  a  function  of  time  in  PI  (5)  and  P2(5)  have  also  been  cor¬ 
rectly  attributed  to  gain  generator  thermal  ring  growth.  Ring  growth  pushes  the  gain  width  outward  to 
better  couple  with  the  mode  width  of  the  resonator  and  away  from  the  obstruction  caused  by  protec¬ 
tive  clippers  at  the  nozzle  exit  plane.8  Likely  both  effects,  a  growth  and  ring  growth,  are  present. 
Lower  J  lasing  lines  such  as  J  =  5  are  more  sensitive  to  this  thermal  growth  because  of  the  sharper 
peaked  spatial  distributions  of  the  low  J  small-signal  gains  upstream  in  the  medium.  However,  for 
HL909B,  Pl(6)  experiences  an  anomalous  “dip”  at  13.5  s  before  recovering.  Thus,  in  fact,  HL909B 
data  are  not  as  convincing  a  validation  as  some  earlier  spectral  data,  such  as  HL909A,8  in  showing 
steadier  power  as  a  approaches  a  steady  value  at  long  operating  times. 

The  third  check  is  on  observed  laser  power  as  a  function  of  time,  P(t),  as  shown  in  Table  2.  In  this 
table,  the  reported  outcoupled  laser  power  from  test  HL909B  is  normalized  to  its  maximum  observed 
value,  between  12  and  14  s.  To  a  zeroth-order-level  approximation,  we  assume  that  the  laser  power 
extracted  must  be  proportional  to  the  F-atom  flow  into  the  resonator,  and  therefore  degree  of  disso¬ 
ciation,  normalized  to  our  computed  quantity,  at  14  s. 

This  normalization  can  be  compared  with  the  actual  similar  trend  in  normalized  data  to  suggest  that 
the  slight  empirical  increase  in  laser  power  as  function  of  time,  P(t),  might  be  due  partly  to  the  slow 
increase  to  steady  state  of  parameter  a.  As  noted  earlier,  ring  growth  also  contributes  to  P(t)  increase, 
and  may  be  the  reason  for  the  larger  part  of  the  P(t)  increase  after  10  s.  Nevertheless  the  time- 
dependent  pressure  curve,  Pc,  plus  slight  temporal  increase  in  reagent  flows  generate  a  performance 
prediction  from  our  model  consistent  with  an  observed  power  trend. 
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Table  2.  Comparison  of  Adjusted  Model  to 
Normalized  HL909B  Outcoupled 
Power8 


Time  (s) 

Normalized  _ 

HL909B 

9 

0.86 

0.71 

10 

0.96 

0.87 

11 

0.92 

0.96 

12 

0.96 

1.00 

13 

0.98 

1.02 

14 

1.00 

1.00 

Outcoupled  powers  are  normalized  to  values  at  t  =  14  s. 


The  general  trends  exhibited  by  the  three  empirical  checks  suggest  that  the  model  is  consistent  with 
observations.  However,  verification  remains  inconclusive  without  more  data  and  experiments. 
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9-  Model  Uncertainty 


The  exponential  structure  of  the  equations  for  F-atom  density,  Eqs.  19-20  combined,  suggest  that  any 
calculation  of  that  parameter  using  empirical  data  will  be  subjected  to  accumulating  errors  in  the  F2  to 
F  density  ratio  in  Eq.  (23).  If  one  supposes  that  the  ratio  of  densities  of  F2  to  F-atoms  can  be 
established  to  the  reasonable  uncertainty  of  ±20%,  then  the  fractional  uncertainty,  Aa/a,  can  be 
shown  to  depend  on  the  value  of  a  as  in  Figure  6.  The  uncertainty  is  computed  assuming  a  root-sum- 
square  fractional  error  dependence  for  independent  error  sources  in  experiments  that  are  not  repeated. 

For  example,  the  uncertainty  in  estimating  a  =  0.9  is  ±2.5%.  On  the  other  hand,  if  a  is  computed  to 
be  0.4,  the  uncertainty  is  ±12%.  The  uncertainty  in  knowing  a  approaches  the  uncertainty  in  estab¬ 
lishing  the  F2/F  density  ratio  for  very  small  values  of  a.  For  any  value  of  a,  the  density  of  F2  is  diffi¬ 
cult  to  establish  even  though  the  total  number  of  F-atoms  can  be  easily  computed  from  thermochemi¬ 
cal  calculations.  There  is  obviously  less  uncertainty  in  this  density  as  a  approaches  unity;  thus,  it  is 
clear  that  the  model  is  most  useful  for  degrees  of  dissociation  approaching  unity. 


Error  in  Degree  of  Dissociation  Computation 
assuming  20%  error  in  density  ratio  [F2/F] 


alpha 


Figure  6.  The  uncertainty  in  computing  degree  of  dissociation,  a. 
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10.  Closure 


In  reality,  accuracy  limitations  of  experimental  laser  data  and  the  chemical  kinetic  database,  espe¬ 
cially  the  high  sensitivity  of  equilibrium  constant  K  to  small  uncertainties,  may  make  it  futile  to  strive 
for  high  accuracy  in  the  density  ratio.  Then  the  utility  of  such  models  for  detailed  design  is  necessar¬ 
ily  limited.  We  contend  that  such  models  are  important,  perhaps  even  invaluable,  tools  for  guidance 
and  insight,  but  may  not  be  very  useful  for  detailed  design  or  planning  except  when  combustor  tem¬ 
peratures  are  set  sufficiently  high,  with  the  degree  of  fluorine  dissociation  nearing  unity. 

The  reader  is  cautioned  that  although  this  model  is  approximately  correct  as  the  combustor 
approaches  “steady  state,”  subsequent  work  generated  a  better  formulation  to  handle  the  large 
unsteadiness  in  flow  parameters  at  the  onset  of  ignition  and  the  approach  to  this  quasi-steady-state 
condition.  It  was  also  later  found  that  the  combustor  stagnation  pressure,  Pc,  could  not  be  success¬ 
fully  modeled  as  a  function  of  time  without  developing  the  concept  of  a  colder  inner  cylindrical  core 
in  this  model.  These  adjustments  are  contained  in  the  latest  and  final  versions  of  the  model. 
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Appendix  A — Mass  Data  for  the  Alpha  Gain  Generator 


The  following  mass  data  was  acquired  from  the  Alpha  Laser  Optimization  program. 


Normal  Ring 

11.6  kg 

End  Ring 

11.0  kg 

Sub  assemblies  per  ring 

0.85  kg 

Domes 

1.15  kg 

preceding  page  blank 
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Appendix  B— Calculation  of  L/e  as  a  Function  of  Run  Time 


In  a  somewhat  circular  argument,  we  assume  L/e  is  nearly  constant,  with  L  =  L0  =  18.2  mm  and  e  =e0 
=  0.6634  mm  defined  by  the  system  geometry.  By  Eq.  (14),  we  can  estimate  AT-GGA.  Using  this 
temperature,  we  can  estimate  the  increase  in  L  (Eq.  13)  and  the  decrease  in  e  (Eq.12)  as  a  function  of 
time. 


t(s) 

po  (psia) 

At  ( 1  s) 

Apo/po=? 

AT-GGA 

e  (mm) 

L  (mm) 

L/e 

4 

30 

0.6638 

18.20 

27.42 

5 

32.8 

4  to  5 

0.0933 

37.0 

0.6621 

18.22 

27.51 

6 

36.1 

5  to  6 

0.1006 

39.9 

0.6603 

18.23 

27.61 

7 

37.7 

6  to  7 

0.0443 

17.6 

0.6595 

18.24 

27.66 

8 

39.8 

7  to  8 

0.0557 

22.1 

0.6584 

18.25 

27.72 

9 

41.7 

8  to  9 

0.0477 

18.9 

0.6576 

18.26 

27.76 

10 

43.5 

9  to  10 

0.0432 

17.1 

0.6568 

18.26 

27.81 

11 

45.5 

10  to  11 

0.0460 

18.2 

0.6559 

18.27 

27.86 

12 

47.1 

11  to  12 

0.0352 

13.9 

0.6553 

18.28 

27.89 

13 

48.9 

12  to  13 

0.0382 

15.2 

0.6546 

18.28 

27.93 

14 

50 

13  to  14 

0.0225 

8.9 

0.6542 

18.29 

27.95 

15 

51 

14  to  15 

0.0200 

7.9 

0.6538 

18.29 

27.98 

From  t  =  4  s  to  t  =  15  s,  L/e  increases  only  around  2%,  thus  validating  our  assumption  that  L/e  is  nearly  constant. 
We  have  noted  through  References  9  and  10  that  estimating  AT-GGA  by  this  method  for  early  times 
may  not  be  exactly  correct  because  time-dependent  changes  in  p0  may  depend  on  flow  variations  as 
well  as  the  thermal  effects  of  this  model. 
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Appendix  C — Estimate  of  Volumetric  F-Atom  Loss  from  Plenum  to  Exit 


We  contended  that  the  volumetric  loss  in  a  one-dimensional  nozzle  flow  from  plenum  to  exit  of  F- 
Atoms  is  around  1%.  The  contour  pattern,  A/A*  for  the  HYWN  nozzle  was  made  available  by  the 
contractor.  This  does  not  include  boundary  layer  effects.  The  table  below  summarizes  our  result  for 
a  high  degree  of  fluorine  dissociation  in  the  plenum.  Parameter  “x”  is  flow  in  the  radial  direction, 
“A/A*”  is  the  nozzle  area  expansion,  “M”  is  the  Mach  number,  “T/T0”  is  the  relative  static 
temperature  at  station  x,  “[F]”  is  the  molar  concentration  of  F-atoms  (Eqs.  19-21),  and  “alpha”  (Eq. 
23)  is  the  degree  of  dissociation.  The  ratio  of  alpha  at  x  =  27.3  mm  and  x  =  0  mm  (plenum)  shows 
less  than  a  1%  loss. 


x(mm) 

A/A* 

M 

T/T0 

[F] 

alpha 

0 

7.41 

0.08 

0.9983 

1.9371  E-06 

0.964 

1.1 

6.48 

0.09 

0.9978 

1.9204E-06 

0.962 

2.2 

4.63 

0.12 

0.9957 

1.9008E-06 

0.960 

3.3 

2.78 

0.21 

0.9878 

1.8652E-06 

0.959 

4.4 

1.57 

0.40 

0.9571 

1.7583E-06 

0.958 

5.7 

1.00 

1.00 

0.7810 

1.222  IE-06 

0.958 

0.0 

1.00 

1.00 

0.7810 

1.2217E-06 

0.958 

1.1 

3.45 

3.04 

0.2790 

1.9502E-07 

0.958 

2.2 

3.98 

3.23 

0.2553 

1.6639E-07 

0.958 

3.4 

6.63 

3.94 

0.1871 

9.5607E-08 

0.958 

4.5 

9.54 

4.48 

0.1507 

6.5001  E-08 

0.958 

5.7 

10.87 

4.69 

0.1396 

5.671  E-08 

0.958 

6.8 

13.25 

5.01 

0.1243 

4.6122E-08 

0.958 

9.1 

16.44 

5.38 

0.1096 

3.6854E-08 

0.958 

11.4 

18.56 

5.60 

0.1021 

3.2509E-08 

0.958 

13.6 

21.21 

5.84 

0.0946 

2.8344E-08 

0.958 

15.9 

23.32 

6.02 

0.0895 

2.5707E-08 

0.958 

18.2 

25.45 

6.19 

0.0851 

2.3487E-08 

0.958 

20.4 

29.16 

6.46 

0.0787 

2.0433E-08 

0.958 

22.7 

30.48 

6.55 

0.0767 

1 .9528E-08 

0.958 

25 

32.61 

6.69 

0.0738 

1.8222E-08 

0.958 

27.3 

34.99 

6.83 

0.0709 

1.6966E-08 

0.958 

In  this  calculation,  we  have  taken  plenum  conditions  to  be  T0  =  1858K,  p0  =  39.8  psia,  y  =  1.56,  and 
W  =  11.1  g-mol-1  as  determined  by  mole  fractions  of  products  (See,  for  example,  Table  1).  Initial 
alpha=  0.964  in  the  plenum. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  specializing  in 
advanced  military  space  systems.  The  Corporation’s  Laboratory  Operations  supports  the  effective  and  timely 
development  and  operation  of  national  security  systems  through  scientific  research  and  the  application  of 
advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and 
its  ability  to  stay  abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and  CCD 
detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid-state 
laser  design,  micro-optics,  optical  communications,  and  fiber-optic  sensors;  atomic  frequency 
standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam 
control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  failure 
analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion;  analysis 
and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle  fluid 
mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and  electric 
propulsion;  environmental  chemistry;  combustion  processes;  space  environment  effects  on 
materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  lubrication  and  surface  phenomena.  Microelectromechanical  systems  (MEMS)  for 
space  applications;  laser  micromachining;  laser-surface  physical  and  chemical  interactions; 
micropropulsion;  micro-  and  nanosatellite  mission  analysis;  intelligent  microinstruments  for 
monitoring  space  and  launch  system  environments. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic-ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric  physics, 
density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmospheric  radiation; 
solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  surveillance,  imaging  and 
remote  sensing;  multispectral  and  hyperspectral  sensor  development;  data  analysis  and 
algorithm  development;  applications  of  multispectral  and  hyperspectral  imagery  to  defense,  civil 
space,  commercial,  and  environmental  missions;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and  magnetosphere;  effects  of 
electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation,  design, 
fabrication  and  test;  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state- specific  chemical  reactions,  and  radiative  signatures  of 
missile  plumes. 


